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Active galaxies are characterized by prominent emission from their nuclei. In the 
unified view of active galaxies, the accretion of material onto a massive compact 
object - now gen erally believed to be a black hole - provides the fundamen- 
tal power source (jReesl. 1984 ). Obscuring material along the line of sight can 



account for t he ob served differences in nuclear emission (iKrolik &i Begelman 



Il988l: Krolikl . 19991 ). which determine the classification of AGN (for example 
as Seyfert 1 or Seyfert 2 galaxies). Although the physical processes of accre 
tion have been confirmed observationally (Greenh ill et al. . 19951: ICallimore et~aL 



19971 ). the structure and extent of the obscuring material have not been deter- 
mined. Here we report observations of powerful hydroxyl (OH) line emissions that 
trace this obscuring material within the circumnuclear environment of the galaxy 
Markarian 231. The hydroxyl (mega)-maser emission shows the characteristics of 
a rotating, dusty, molecular torus (or thick disk) located between 30 and 100 pc 
from the central engine. We now have a clear view of the physical conditions, the 
kinematics and the spatial structure of this material on intermediate size scales, 
confirming the main tenets of unification models. 

Observations of powerful molecular water-vapour and hydroxyl megamaser 
emissions , which are a million times more luminous than their galactic coun- 
terparts (Baan, 1985[ l. have provided unique information about distinct re- 
gions within galactic nuclei. Water-vapour masers can map out the central 
few parsecs of an AGN, tracing the rotation of nuclear a ccretion disks a nd the 
state of molecular material surrounding the radio jets ( Malonev . 2002| ). The 
inferred geometry and thickness of these disks, however, cannot account for 
the observed obscuration signature and the radiation patterns seen in active 
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galactic nuclei (AGN) (|Malonevl . liooi iHerrnstein et al.l . Il998| ) . On the other 
hand, t he hydroxyl masers in galactic nuclei are spread across hundreds o f 
parsecs (jPiamond et all Il999l : [Lonsdale et aP . Il998': 'Pihlstrom e t allboniL 
Their tight correlation with the infrared emission fj.Baan . .1989.) indicates a 
strong association wi th the dusty obscuring material along t he line of sight 
towards the nucleus ( Krolik &: Beeelman . 19881 : Krolik . 1999| ). To study the 
distribution of the hydroxyl megamaser emission with respect to the nuclear 
engine at a resolution of tens of parsecs in Mrk 231, we performed very-long- 
baseline interferometry (VLBI) radio observations with the European VLBI 
Network (EVN). Mrk 231 is the most luminous infrared galaxy in the local 
universe (distance of 172 Mpc), and is undergoing a merger as indicated by 
the high dust content and the prominent tidal tails obs erved at optical wa ve- 
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lengths. The nuclear power source has been classified ()Baan et al 
a Seyfert 1, which in the unification scheme would suggest that the observer 
has a direct, and almost unobscured, line-of-sight view into the broad-line 
emission region at parsec scales. The location of the nuclear engine on simi- 
lar scale sizes is traced by the radio outflow (lUlvestad et al. . 1999b[ ). which 
extends to almost a hundred parsecs (|Ulvestad et al ■L ll999al) . 

The EVN data reveal an OH emission structure with an east-west extent of 
130 pc that is straddling a north-south elongated radio continuum emission 
(Fig.P). This c ontinuum emiss i on has been further resolved into a lobe-core- 



lobe structure (|Ulvestad et al.L |l999a|) of 80 pc at position angle PA = 8° 



which i s found to be rn i saligne d by 57° with the core radio structure on parsec 
scales f Ulvestad et al. . 1999bl ). The emission spectrum integrated across the 
whole OH emitting region displays the two main hydroxyl lines at 1667 and 
1665 MHz showing a line-width of around 200 kms~^ (full-width at half- 
maximum, FWHM) at a centroid velocity of 126 10 kms~VFig.ll). Th e EVN 
data and also the intermediate-scale MERLIN ( Richards et al.L 2000| ) obser- 
vations account for only half of the integrated line emission seen by the West- 
erbork Synthesis Ra dio Tel e scope (WSRT; see Figure 2) and previous single- 
dish measurements f Baan . 1985| ). This missing emission component must 
therefore originate within the resolution gap between the EVN and WSRT 
observations, in an extended region between 190 pc and 11 kpc of the nucleus. 
This component is p robably associated with the nu clear continuum ernission 
(iTavlor et al .'.'l999l) and t he CO(l-O) line emission (jBrvant fc Scovill3 . ll996l : 



Downes fc Solomon! 1998[ ). Figure 3 (top panel) presents the spatial velocity 



distribution, which is similar for each of the two hydroxyl lines and which has 
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an overall velocity shift of about 160 kms~^ from the northwest towards the 
southeast. The velocity width of the emission lines has an average value of 
about 70 kms~^ (FWHM) across the northern part of the radio continuum 
emission. On the other hand, there is also a distinct OH emission region at 
the northwestern edge of the continuum emission that has a velocity similar 
to that seen towards the centre of the radio core, and which has a distinctly 
higher velocity width of 85 kms~^ (Figure 3; middle panel). 

The nature of the OH emission observed in the nucleus of Mrk 231 adds 
a new component to the hierarchical structure of this system, and comple- 
ments the information obtained with other tracers at lower resolution. We 
find that the OH velocity field with its symmetry axis at PA = 34°does not 
agree with that of the kiloparsec- scale CO(l-O) structure ( Brvant fc Scovill3 . 
19961 : Downes fc Solomonl . 1998| ) having a PA = 0°. Neither does it agree 
with the large- scale radio outflow (PA = 8°), or the more compact C0(2- 
1) distribution ( Downes fc Solomon . 1998[ ) extending up to 850 pc at PA = 
18°. However, the OH velocity field is in agreement with that of the line- 
of-sight HI absorption (ICarilh et al.L Il998t ) at PA = 27°. The HI absorption 
is detected against a diffuse radio-emitting halo with an extent of about 
350 pc and an orientation of PA = 22°. This diffu se radio halo, which is pos- 
sibly associated with a circumnuclear starburst (jTavlor et al. ,119991), could 
then also trace th e dense and opaque interstellar medium of a disk or torus 
1998) inclined at 56°. The integrated HI velocity dispersion 
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of 193 ± 25 kms~^ agrees well with that of the hydroxyl 
lines (Fig. |2)). This suggests that the OH emission indeed originates in the 
central regions of a circumnuclear rotating disk or torus having an outer 
structure that is traced by the C0(2-l) emission. 

The energy source responsible for the excit ation of the OH molecules is the 
extreme far-infrared radiation ( Baan . 1989| ) field in Mrk 231. The observed 
far-infrared spectral temperature of around 43 K provides the ideal pumping 
conditi ons for the OH clouds, which could ran ge in size from 10 pc up to 



100 pc ( Henkel et al. . 1987 : Randell et al. . 1995| ). The ratio of the intensity of 



the two maser lines varies rather smoothly across the entire emission structure 
and has an average value of 1.8, as expected for local temperature equilibrium 
(LTE). This would indicate that the observed emission lines are optically 
thin and are mostly unsaturated masers. Assuming that the diffuse radio 
continuum, which serves as a background for the nuclear HI absorption, also 
serves as a background for the observed OH masers, this emission could 
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be explained in terms of exponential amplif ication with a maximal gain o f 
about 2.2 within classical OH maser-models ( Baan et al. . 19821 : Baan . 1985| ). 
However, extreme values for the line ratio have been found in the discrepant 
northwestern region discussed earlier in this text. This region must have a 
completely different pumping environment than the rest of the structure. The 
unusual blue-shifted velocity, the slightly higher velocity dispersion, and the 
line ratios o f this discrepant region are all consistent with the theoretically 
predi ctions ( Baan . 19851 : Randell et al. . 1995t Wada fc Norman . 20021 : Baan . 
Si) for an interaction region between the radio outflow and the molecular 
environment. 

The molecular environment described above for the nuclear region of Mrk 231 
is very similar to the dust structure s found in active galax i es suc h as NGC 4261 
using the Hubble Space Telescope ( van Langevelde et al. . 2000l ). We have ap- 
plied an axisymmetric torus model with a circular cross-section to the EVN 
data, in order to determine the structural boundaries in the nuclear region. 
The compact molecular structure surrounding the nucleus is misaligned with 
the large-scale molecular disk of the galaxy by 34°. The radius of the in- 
ner cavity of the torus would be about 30 pc, as suggested by the location 
of the northwestern interaction region. The outer edge of the torus spans 
a diameter of 200 pc, which is ba sed on the 4 mJy contours of the diffuse 



radio emission dCari 



li et al.L 119981 ) and on far-infrared blackbody estimates 



( Downes fc Solomon . 11998^ This region also serves as a radio background 



for the OH amplification process. These parameters result in a torus struc- 
ture with a radius of 65 pc and a thickness of 70 pc, which accounts for 
an obscuration angle of 60° centered on the plan e of the torus and wh ich is 



representative of values found for other galaxies (jSchmitt et al.Ll2001[ ). The 



orientation and the shape of the extended radio contours and the symme- 
try axis of the velocity field indicate that the torus is tilted upward from 
the line of sight by 56°, and is rotated anticlockwise by 35°. This inferred 
model for the torus has been displayed as a wire-structure representation 
and is presented in Figure 3 (bottom panel). The nuclear radiation cones 
with opening angle s of 60° , assumed to be similar to those seen in M 87 



( Junor et al. . 19991 ). represent the directions with an unobscured view of the 



nucleus. The radio outfiow follows a twisted path within these cones, starting 
at the nucleus (at PA = 65°) and ending up as a d ouble source aligned w ith 
the symmetry axis of the outer disk (at PA = 8°) (|Ulvestad et al. . 1999bl a). 
As a result, this parameterization indicates that the nuclear accretion disk. 
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that coUimates the jet outflow inside the central cavity, is actually misaligned 
by 30° with the plane of the molecular torus. 

In reahty, the proposed torus structure does not need to have a circular 
cross-section. In addition, it will most probably be wrapped inside a cocoon- 
like surface region with a higher temperature, and will have the same outer 
extent of 460 pc as the diffuse radio emission and the HI absorbing gas 



)U p 

( Carilli et all EqQS). However, the representation of the inner part of this 



rotating dusty and molecular torus in Mrk 231 is consistei it with ( all) cu rrent 
models of galactic nuclei and theoretical investigations ( KrolikL [l999[ ). and 



supports the uniflcation schemes for AGN. Whether a different orientation 
of the torus and the nuclear accretion disk represents a special case or is 
a general characteristic of active nuclei needs to be investigated with VLBI 
observations of other megamaser galaxies hosting a similar nuclear power 
source. 
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Figure 1: Hydroxyl and radio— continuum emission in the nucleus of 
Mrk 231. 

The integrated OH-hne emission in pseudo-color (in mJy per beam) is su- 
perimposed on the nuclear continuum emission in contours (peak 39 mJy 
per beam). The contour levels are a geometric progression in factors of the 
square root of 2 starting at 4 mJy per beam. These EVN observations of 
Mrk 231 were made on September 1999 for 12 hrs in dual-polarization phase- 
referencing mode at 1599 MHz and with 256 spectral channels across a total 
bandwidth of 8 MHz. The synthesized beam is 39 mas x 28 mas, where 
1 mas corresponds to a size of 0.83 pc at the distance of Mrk 231 (172 Mpc 
by assuming qo=0.5 kms~^ and Ho=75 kms~^ Mpc~^). The data has been 
correlated at the EVN correlator at the Joint Institute for VLBl in Europe 
(JIVE), and later transferred to the AlPS packages for calibration of the 
complex visibilities. The integrated line emission image has been obtained 
by averaging over the velocity range of the OH lines, after subtraction of the 
radio continuum, obtained by averaging the channels without line emission. 
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Figure 2: Hydroxyl emission spectra of Mrk 231 traced at different 
scale sizes. 

The top panel shows the WSRT spectrum taken at 14 arcsec resolution, the 
middle panel displays the spectrum for the EVN observations at 39 mas 
resolution, and the bottom panel displays the difference spectrum of both 
observations. The emission line features show both OH main lines at 1667 
and 1665 MHz that are offset by 365 kms^^; the velocity scale refers to 
the 1667 MHz emission line by assuming a heliocentr ic velocity of Vrpntpr 
= 12 650 kms~^ based on the CO(l-O) centre velocity (jPownes fc Solomon . 
IQQSi l. The spectral resolution is 18 kms~^ per channel. The OH line width 



in the EVN data is 214 ± 11 kms"^ (FWHM), which is 70 kms"^ less 
than for the WSRT data. The residual emission line in the bottom frame 
counts for 44 % of the OH emission in Mrk 231 and is systematically red- 
shifted by around 45 kms~^. Similar flux discrepancies, but with no dis- 
tinct velocity o ffset, have been found for other OH-MM galaxies studied at 
high resolution ( Diamond et al.> . 1999:. Lonsdale et al. . 19981 : IpThlstrom et al. 



20011). Similar velo c ity off sets have also been found for the CO ( 2-1) emission 
flSrvant &: Scoviil3 . 1996t ) and the HI absorption in Mrk 231 (|Carilli et al 



1998t) . 




Figure 3: The circum— nuclear kinematics in Mrk 231. 

Top panel - The velocity field of the OH 1667 MHz emission in the EVN 
data is shown in pseudo-colors. The spatial resolution is the same as for 
Fig. n and the velocity scale is the same as in Fig. El The overall velocity 
pattern reveals an east-west velocity gradient of 1.41 kms~^ per parsec over 
the entire OH emission such that the northwestern edge of the torus moves 
towards the observer. 

Middle panel - The line width of the 1667 MHz emission lines (in kms~^) is 
superimposed on the nuclear radio continuum emission in contours displayed 
in a similar way as in Fig. Q The diagram shows that the velocity dispersion 
varies significantly across the emission region. The velocity dispersion struc- 
ture shows several distinct regions but with the highest values to be found 
at the western edge of the radio continuum source. 

Bottom panel - The inferred model of the nuclear torus in Mrk 231 is dis- 
played as a wire diagram with symmetric ionization cones. This model takes 
into account all large-scale characteristics of the nuclear radio emission and 
the OH emission. The symmetry axis of the torus has been inclined upwards 
by 56°, and then rotated anticlockwise by 35°. The molecular material moves 
from top-right to bottom-left (northwest to southeast). The virial estimates 
of the central mass concentration is 7.2 ± 3.8 xlO^ M^. 



